Multielectron ionization of atoms in strong fields: 
Classical analysis in symmetric subspace 
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We consider the final stage of multiple ionization of atoms in a strong linearly polarized laser 
field within a classical model. We propose that non-sequential multiple ionization is dominated by 
symmetric escape from a highly excited intermediate complex. For a configuration of TV electrons 
with Cnv symmetry in a plane perpendicular to the electric field we analyze the classical motion 
in phase space and discuss the final momentum distribution parallel and perpendicular to the po- 
larization axis. The results are in good agreement with recent experiment of multiple ionization in 
Ne. 
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Double (or multiple) ionization of atoms is a funda- 
mental process to our understanding of many electron 
dynamics in external fields. The surprisingly large yields 
of multiply charged ions reported in the first experiments 
with intense laser fields |Q clearly show that sequen- 
tial ionization is not the leading process and that the 
electron-electron correlated dynamics has to be taken 
into account. One possibility is that the electron is 
shaken-off due to a non-adiabatic change in potential 
during the ionization of the first electron a process 
that accounts for double ionization at high photon en- 
ergies (above 1 keV) j|. For multiphoton ionization a 
two step process is more likely: one electron is ionized 
first, accelerated by the field and driven back to the atom 
where it ionizes the second electron in a rescattering col- 
lision |^-|lC| ] . The applicability of tunneling expressions 
as in EDO] then is due to the first ionization taking place 
through a tunneling process in a quasi-static electric field. 
Such a sequence of events is supported also by the nu- 
merical calculations of Becker and Faisal . 

The correlation between the electrons was convincingly 
demonstrated in a series of recent experiments. Measure- 
ments of the distributions of ion and electron momenta in 
double and triple ionization |p^L5| clearly show a pref- 
erence for electrons escaping towards the same side of the 
atomic core along the field polarization axis 

0- On first 

sight this seems to be very different from the symmetric 
escape of electrons in double ionization in the absence of a 
field, where according to Wanniers analysis the electrons 
escape in opposite directions Jl6],[l7|]. A closer analysis 
of Wanniers arguments shows, however, that they can be 
applied to this situation as well p8| . It is our aim here 
to generalize this model to the correlated escape of three 
or more electrons. 

The essential elements of Wanniers analysis are a di- 
vision of the process into two steps, the formation of a 
highly excited complex of two electrons close to the core 
and the two electron escape from this complex. At the 
threshold for ionization, energy is scarce. Therefore, in 
the dominant channel, mutual repulsion should be mini- 



mal and the energy should be equally shared by the out- 
going electrons. Any deviation from the symmetric ar- 
rangement in phase space or in energy would be amplified 
and lead to single rather than multiple ionization. These 
requirements become less stringent the higher the energy 
above threshold. 

One can argue that a similar division of the ioniza- 
tion process is possible in electric fields |lq] . In partic- 
ular, ionization also proceeds through the formation of 
an intermediate highly excited complex, created during 
the rescattering event. But the next step is different, 
the configuration that dominates the ionization channel 
is no longer the Wannier state and the threshold energy 
is modified as well. For instance, double ionization is 
observed even though the estimated energy transferred 
by returning electron is too small for immediate double 
ionization As discussed in |L8| this is possible if 

during the collision the external field is not zero, for then 
a Stark saddle opens through which the electrons can 
escape. This saddle breaks the symmetry and focuses 
the electrons in the direction of the electric field. It also 
forces the electrons to be close to a symmetry subspace, 
since by mutual repulsion the electron that reaches the 
saddle first can push the other back to the nucleus. This 
would then result in either single ionization or in another 
rescattering event, but not in double ionization. Classi- 
cal trajectory calculations within this symmetry subspace 
are able to reproduce the main features of the experimen- 
tally observed ion momentum distributions. 

In the present letter we present a generalization of our 
model 1 18 1 to multielectron ionization. The key assump- 
tion is that the process is dominated by a symmetric 
configuration of the electrons with respect to the field po- 
larization axis as suggested by the experiments [|12|-|l5[. 
Specifically, we assume that all electrons move in a plane 
perpendicular to the field and that they obey a Cjv« sym- 
metry, which generalizes the Ci v symmetry of the previ- 
ously analyzed case. The reflection symmetry limits the 
momenta to be parallel to the symmetry planes and thus 
confines the motion to a dynamically allowed subspace in 
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the high-dimensional TV-body phase space. I.e. if in the 
full phase space of the TV-body problem initial conditions 
are prepared in this subspace they will never leave it. 

We take the electric field directed along the z-axis, 
and the positions of the iV-electrons as Zi = z, pi = p 
and ifi — 2ni/N, where (pi,ifi,Zi) are cylindrical coor- 
dinates. The momenta of the electrons are all identical, 
P P ,i = P P , Pz.i = Pz and p v ,i = 0. For this geometry 
the classical Hamiltonian for N electrons, for zero total 
angular momentum along the field axis, in atomic units, 
with infinitely heavy nucleus and in dipole approxima- 
tion, reads 

H(p p ,p z ,p,z,t) = N^^ + V(p,z,t), (1) 



with potential energy 

iV 2 N(N-l) 



V 



Vp 2 + z 2 4psm(n/N) 



+ NzFf(t) cos(tot + , 



and the pulse shape 

f{t) =sin 2 (7rf/T rf ). 



(2) 



(3) 



F, to and </> stand for peak amplitude, frequency and 
phase of the external field respectively, while Td is the 
pulse duration. 

Experiments show that the multiple ionization is possi- 
ble even when the energy transferred by the rescattered 
electron is less than that needed to lift the system to 
the many electron continuum. We therefore assume that 
the complex from which multiple ionization starts is de- 
scribed by Hamiltonian ([!]) with initially negative energy 
E. Moreover, since the electronic motion close to the nu- 
cleus is much faster than the changes in the electric field 
the process of electron escape can be discussed adiabat- 
ically for fixed external field. This allows us to analyze 
the potential energy for a fixed time. 

Before we proceed with the analysis we take advantage 
of the scaling symmetry of the classical Hamiltonian (fil) 
and eliminate one parameter. If the variables are rescaled 
according to 



H -> VFH, 
p, z -» p/VF, z/Vf, 
p p , Pz -> F l /%, FV*p z 
u, t -» F 3 / 4 u, F-^H 



(4) 



the dynamics becomes independent of the peak value of 
the field amplitude, i.e. we obtain the system described 
by the Hamiltonian ([!]) with F = 1. 

Equipotential curves of Eq. (||) in the scaled variables 
(@), for N = 3 and N = 6 are shown in Fig. [|. The 
saddles are located along the lines z — r s cos 9 S and 
p = r s sin 9 S with 9 S = 9 or 9 S = it — 9 where 



9 = arctan 



47Vsin(7r/7V) 
N - 1 



2/3 



-1/2 



and 



N 



N - 1 



2/3 



(5) 



(6) 



s ft V \4:Nsm(n/N) 
with ft = \f(t) coa(ujt + <j>)\. The energy of the saddle is 



V s = -y/4N3f t 



1 - 



N — 1 



47Vsin(7r/jV) 



2/3' 



3/4 



(7) 



During a field cycle the saddle moves in from infinity 
along the line 9 S — 9, back out to infinity and then in 
and out again along the line 9 S = n — 9. The angle 9 
increases with the number of electrons N, while the en- 
ergy of the saddle changes non-monotonically as shown in 
Fig. ^. For N > 14 the saddle disappears. For this many 
electrons the repulsion between the electrons is stronger 
than the attraction to the nucleus, see Eq. (0) for z = 0. 
The situation with N > 14 is rather an extreme case 
from the experimental point of view but it reflects the 
fact that a bound state for a neutral atom with N elec- 
trons symmetrically distributed in a plane can not be 
formed if N > 14. 

The saddle which is opened by the external field al- 
lows electrons to escape even if the total energy E of the 
highly excited complex is negative. It is enough, classi- 
cally, that E is greater than the energy of the saddle. The 
saddle thus provides a kind of transition state jl9 20 for 
the ionization process: once the electrons cross it, they 
are accelerated by the field and pulled further away, mak- 
ing a return rather unlikely. Moreover, they can acquire 
the missing energy so that the electrons can escape even 
when the field vanishes. The field thus plays a double role 
in determining a threshold for this process: during the 
first stages of the rescattering process it provides the en- 
ergy for the collision complex and during the final stages 
it opens the path for multiple escape. 

In the experiment p]] of triple (also single and dou- 
ble) ionization of Ne by ultrashort (30 fs) laser pulses at 
795 nm and at intensities 1.5 PW/cm 2 the distributions 
of momenta of ions parallel and perpendicular to the po- 
larization axis have been measured. In the limit of negli- 
gibly small momentum transfer by the absorbed photons, 
the ion momentum pi on reflects the sum of the momenta 
of the emitted electrons, pi on = —J2iPi P^Jl^] - To cal- 
culate the experimental distributions corresponding to 
the triple ionization of Ne we have performed classical 
trajectory simulations for N — 3. Note that the correct 
equations of motion for a single electron follow from the 
Hamiltonian (]l|) after division by N. 

The Hamiltonian ([!]) with N = 3 is the full Hamilto- 
nian for a Li with its three electrons. In order to relate it 
to triple ionization from atoms with more electrons the 
zero point in energy has to be shifted to the threshold for 
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the three electron continuum, and all energies have to be 
taken relative to this level. Moreover, interactions with 
the electrons in the core are neglected. Specifically, for 
the modeling of the experiments |15f| on triple ionization 
in Ne we thus assume that in a rescattering process the 
energy transfer is less than the threshold for triple ion- 
ization (about 4.6 a.u.). The precise value of the initial 
energy E depends on the details of the rescattering pro- 
cess and can thus not be determined within our model: it 
is a free parameter. The rescattering event is most likely 
to take place when the field amplitude is high, so we have 
started the simulation close to the top of the pulse, i.e. 
for to — 0.33^, where T4 — 2412 a.u. corresponds to 
the experimental pulse duration . The corresponding 
amplitude and frequency of the field are F = 0.207 a.u. 
and uj = 0.057 a.u. For such field parameters the sad- 
dle energy equals —3.5 a.u. We have performed numer- 
ical simulations for several initial energies. The distri- 
butions of the electron momenta for E = — 1 a.u. are 
shown in Fig. || The distribution for parallel momenta 
can be compared with the corresponding experimental 
distribution of ions momenta The main features 

of the distribution, its width and characteristic double 
hump structure, are well reproduced in our calculation. 
The minimum in the distribution is, however, less pro- 
nounced than that observed in the experiment. Overall, 
we take this good agreement as strong indication that 
triple ionization can only occur in the neighborhood of 
the symmetric process discussed here. 

For initial energies below about —3 a.u. the double 
hump structure disappears and only a single maximum 
remains. This corresponds to experiments with weaker 
fields which transfers less energy to the rescattered elec- 
tron. Then the electrons crossing the saddle have smaller 
kinetic energy and the interaction with the remainder of 
the pulse blurs the distribution. Such a change in the 
character of the distribution with decrease of field inten- 
sity has also been observed in the experiment of double 
ionization of He atoms |Q . 

The distribution for perpendicular ion momenta can 
not be calculated in our model since this component van- 
ishes exactly by symmetry. But we can calculate the 
transverse momentum of an individual escaping electron, 
as shown in Fig. ||. The absence of events around zero re- 
flects simply the effect of the repulsion between electrons. 
Measurements of this quantity should provide further test 
of this model. 

In conclusion we have shown that the double humped 
structure in ion momenta distribution observed in the 



We would like to thank Harald Giessen for stimulat- 
ing our interest in this problem and for discussions of 
the experiments. Financial support by the Alexander 
von Humboldt Foundation and by KBN under project 
2P302B00915 are gratefully acknowledged. 



experiment of triple ionization of Ne 15 can be calcu- 
lated within a classical simulation assuming symmetri- 
cally escaping electrons. This multi-ionization path can 
be extended to include up to 13 electrons. The successful 
comparison with the experimental data supports the idea 
that the dominant contribution to ionization comes from 
initial conditions in the highly excited collision complex 
that are asymptotic to this symmetry plane, very much 
as in Wanniers analysis |n||l7|,^l) . 
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FIG. 1. Equipotential curves of the potential energy (|2j), 
in the scaled variables (^J) , for two different values of the num- 
ber of electrons: N = 3 [panel (a)] and iV = 6 [panel (b)]. 
The potentials are plotted for fixed time corresponding to the 
maximal field amplitude. The saddles move along the dashed 
lines when the electric field points in the positive ^-direction 
and along the second obtained by reflections on z = during 
the other half of the field cycle. Note that the angle between 
the z-axis and the dashed line increases with the number of 
electrons. 
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FIG. 3. Final distribution of electrons momenta parallel 
[panel (a)] and perpendicular [panel (b)] to the field polar- 
ization axis for E = — 1 a.u. Panel (a) corresponds to the 
experimental distribution of ions momenta measured in the 
triple ionization of Ne |15| . Note that, in panel (b), the dis- 
tribution rapidly falls to zero around p p w reflecting the 
effect of repulsion between the electrons. The distributions 
are based on an ensemble of 2 ■ 10 5 trajectories. 
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